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ABSTRACT 

We have performed a multi-dimensional correlation analysis on the observed properties of a statistically 
complete core-selected sample of compact radio-loud active galactic nuclei, based on data from the VLBI 
Space Observing Programme (Paper I) and previously published studies. Our sample is drawn from the 
well-studied Pearson-Readhead (PR) survey, and is ideally suited for investigating the general effects 
of relativistic beaming in compact radio sources. In addition to confirming many previously known 
correlations, we have discovered several new trends that lend additional support to the beaming model. 
These trends suggest that the most highly beamed sources in core-selected samples tend to have a) high 
optical polarizations; b) large pc/kpe-scale jet misalignments; c) prominent VLBI core components; d) 
one-sided, core, or halo radio morphology on kiloparsec scales; e) narrow emission line equivalent widths; 
and f) a strong tendency for intraday variability at radio wavelengths. We have used higher resolution 
space and ground-based VLBI maps to confirm the bi-modality of the jet misalignment distribution for 
the PR survey, and find that the sources with aligned parsec- and kiloparsec-scale jets generally have 
arcsecond-scale radio emission on both sides of the core. The aligned sources also have broader emission 
line widths. We find evidence that the BL Lacertae objects in the PR survey are all highly beamed, 
and have very similar properties to the high-optically polarized quasars, with the exception of smaller 
redshifts. A cluster analysis on our data shows that after partialing out the effects of redshift, the 
luminosities of our sample objects in various wave bands are generally well-correlated with each other, 
but not with other source properties. 

Subject headings: galaxies : jets — galaxies : active — quasars : general — radio galaxies : continuum 
— methods : statistical 



1. INTRODUCTION 

One of the major obstacles to a better understanding 
of active galactic nuclei (AGNs) is the disentanglement of 
orientation and relativistic beaming effects on their ob- 
served properties. The use of high-resolution very long 
baseline inteferometric (VLBI) techniques has led to enor- 
mous progress in this regard, but the determination of 
viewing angles and intrinsic speeds of relativistic jets in 
compact extragalactic radio sources remains exceedingly 
difficult. Considerable effort has gone into the study of 
individual objects (e.g., 3C 279: Wehrle et al. 2001; 3C 
345: Lobanov & Zensus 1999), but it is unclear how repre- 
sentative these sources are of the overall radio-loud AGN 
population. As a result, the degree to which relativistic 
beaming biases our view of compact radio sources is not 
well known. 

An effective method of addressing these issues is through 
studies of large well-selected samples of objects, of which 
several are currently underway (e.g., Taylor et al. 1996; 
Kellermann et al. 1998; Fomalont et al. 2000). In a com- 
panion paper (Lister et al. 2001; Paper I), we describe 
one such sample drawn from the Pearson-Readhead (PR) 
survey for which we have obtained high-quality 5 GHz 
space- VLBI images with the VLBI Space Observing Pro- 
gramme (VSOP). The PR survey is arguably one of the 
best-studied radio samples in existence, with an enormous 
amount of supporting data available at a variety of wave- 



lengths. As such, it is ideally suited for studying the effects 
of relativistic beaming on a wide variety of source prop- 
erties, and for testing general predictions of the beaming 
model. 

Here we carry out a multi-dimensional correlation anal- 
ysis of the parsec- and kiloparsec-scale properties of a flat- 
spectrum, core-selected subset of the Pearson-Readhead 
survey. We identify several new trends involving radio 
morphology and intraday variability which suggest that 
these properties are highly influenced by relativistic beam- 
ing. We also confirm previous results which show that 
the core dominance, optical polarization, and emission line 
equivalent width are useful statistical beaming indicators. 
In §2 we describe our sample selection criteria, followed 
by a description of the observed source properties in §3. 
We describe our statistical tests in §4. In §5.1 we use 
cluster analysis to investigate natural groupings of source 
properties, and discuss the role of emission line strength 
as a beaming indicator in §5.2. In Sections 5.3 and 5.4 
we describe Monte Carlo beaming simulations and possi- 
ble ways to identify highly beamed sources. Finally in §6 
we summarize our results and make suggestions for future 
work. 

2. SAMPLE DESCRIPTION 

In Paper I we described a subsample of 31 compact ob- 
jects from the Pearson-Readhead survey that were cho- 
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and define the spectral index according to S v oc v a . 

3.1. General 'properties 

In Table 2 we list various general properties of our 
sources. In columns 1 and 2 we give the IAU and com- 
mon name of each source, and indicate whether it is a 
member of the third EGRET catalog of gamma-ray loud 
sources (Hartman et al. 1999). In column 3 we give the op- 
tical classification as a radio galaxy (RG), BL Lac object 
(BL), or quasar (Q). If a particular quasar has a published 
optical polarization greater than 3% at any epoch we clas- 
sify it as a high-optically polarized quasar (HPQ). For the 
remaining quasars, we classify those with three or more po- 
larization measurements with m opt < 3% as low-optically 
polarized radio quasars (LPRQ). Our BL Lac classifica- 
tions are based on the catalog of Padovani & Giommi 
(1995). 

We have used previously published arcsecond-scale radio 
images to classify our sources into four basic morphological 
types depending on whether they possess distinct hotspots 
on one or both sides of the core (column 4). Those sources 
with no hotspots are termed either compact or halo, with 
the latter displaying weak diffuse structure surrounding 
the core component. 

3.2. Variability properties 

Allcr et al. (1992) have monitored the long-term vari- 
ability properties of the PR sample, and defined a peak- 
to-trough variability index V: 



sen as suitable targets for space- VLBI. These objects all 
have flux densities > 0.4 Jy on the longest Earth baselines 
(<~ 10, 000 km) at 5 GHz. This sample (hereafter referred 
to as VSOP-PR) is well-suited for ensuring likely fringe- 
detection on space baselines, but is not statistically com- 
plete from the standpoint of relativistic beaming. Since an 
accurate relation between intrinsic core and lobe emission 
has not yet been established, current beaming simulations 
can only make predictions for the properties of AGN core 
components. Therefore, one needs a sample that is se- 
lected on the basis of core flux only. The simplest way to 
obtain such a sample is to use a high selection frequency 
(e.g. ;> 15 GHz), at which the steep-spectrum lobe emis- 
sion makes a negligible contribution to the total flux den- 
sity. However, currently available all-sky radio surveys 
only range up to 5 GHz. As a result, nearly all large radio 
surveys (such as the Pearson-Readhead sample) tend to 
contain a mixture of core- and lobe-dominated objects. 

A common method of excluding the lobe-dominated ob- 
jects from radio surveys is to adopt a simple two-point 
spectral-index criterion, such that all objects with spectra 
steeper than a certain value (usually a- < — 0.5 ; S oc v a ) 
are excluded. This method has some drawbacks, as the 
spectral indices may not always be based on simultane- 
ous flux density measurements, and can be affected by 
source variability. Also, some sources whose spectra peak 
around a few GHz (the so-called Gigahertz Peaked Spec- 
trum, or GPS sources) may have flat two-point spectral 
indices, and will be included in the final sample. Recent 
evidence has indicated that these GPS sources, along with 
a similar class, the compact steep-spectrum sources, are 
much younger and more luminous than typical compact 
AGNs, and have weak parsec-scale core components that 
are not likely to be highly beamed (Taylor et al. 1996). 
On this basis they should be excluded from core-selected 
samples. 

In consideration of these issues, we have constructed a 
core-selected subset of the Pearson-Readhead survey by 
using time-averaged 5-15 GHz spectral indexes from the 
University of Michigan monitoring program, as tabulated 
by Aller et al. (1992). These were computed from monthly 
averages of measurements taken over several years, and are 
less susceptible to the variability errors described above. 
By adopting a cutoff of 0:5-15 > —0.4 we effectively elimi- 
nated all of the GPS and compact steep spectrum sources 
from the PR sample. Our final flat-spectrum Pearson- 
Readhead sample (hereafter referred to as the FS-PR) con- 
tains 32 objects, 28 of which are members of the VSOP- 
PR sample. The four additional objects are the quasars 
0723+679, 0850+581, 0954+556, and 2351+456, and are 
targets of upcoming space- VLBI observations with VSOP. 

3. OBSERVED SOURCE PROPERTIES 

We have gathered a vast amount of data on the FS-PR 
and VSOP-PR samples from our space VLBI data pre- 
sented in Paper I and previously published studies. A 
summary of these quantities is given in Table 1, and a 
detailed explanation follows in §3.1. Throughout this pa- 
per we use a standard Freidmann cosmology with H a = 
65 km s _1 Mpc -1 , q a = 0.1, and zero cosmological con- 
stant. We give all position angles in degrees east of north, 

4 The Very Large Array (VLA) is maintained by The National Radio Astronomy Observatory, which is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 
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where S max and S ra j„ are the highest and lowest flux 
measurements, and the cr's are the associated measure- 
ment errors. We will refer to their published indices at 4.8 
and 14.5 GHz in our study as V5 and V15 respectively. 

A large number of Pearson-Readhead sources have also 
been monitored for intraday variability (IDV) with the 
100m Effclsberg telescope and/or the VLA 4 (Quirrenbach 
et al. 2000). These authors have classified sources into 
three distinct IDV classes on the basis of their structure 
functions. Type II objects exhibit significant up-and-down 
variations on hourly timescales, while type I objects show 
monotonic rises or decays in flux on timescales less than a 
day. Type objects show no apparent IDV activity. Sta- 
tistical studies have shown that a particular source may 
have different IDV classifications at different times. For 
those PR sources that have been monitored on one or 
more occasions for IDV, we list in column 5 of Table 2 the 
highest IDV classification from Krichbaum et al. (1992), 
Quirrenbach et al. (1992), and Quirrenbach et al. (2000). 

3.3. Optical spectral properties 

Lawrence et al. (1996) have obtained high-quality op- 
tical spectra of all the Pearson-Readhead survey objects 
with the 5m Palomar telescope. In column 6 we list their 
measured redshifts, and in column 7 we list the rest-frame 
equivalent width of the widest permitted emission line in 
A, corrected to the source frame. 
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3.4. Luminosity and core properties 

In column 8 of Table 2 we list the 5 GHz luminosity of 
the parsec-scale core component, measured in the source 
rest frame assuming a core = 0. The core flux density and 
corresponding rest frame brightness temperature are given 
in Paper I. 

In column 9 we list the total luminosity of the source at 
5 GHz, also in the source frame. We determined the to- 
tal flux density and 5-15 GHz spectral index at the VSOP 
observation epoch (where applicable) by interpolating the 
University of Michigan radio light curve 5 for each source. 
For those sources where 15 GHz data weren't available, we 
used the time-averaged spectral index listed in Aller et al. 
(1992). 

The relative prominence of the unresolved VLBI core 
with respect to extended emission in an AGN is a useful 
indicator of relativistic beaming, since the latter emission 
is thought to be largely isotropic and not highly beamed. 
There is no standard way of defining this parameter in the 
literature, with different authors adopting various ratios of 
core-to-total flux, core-to-cxtended flux, or total parsec- 
scale to total (single-dish) flux. In this paper we adopt 
a parsec-scale core dominance (R pc ), which we define as 
L core /(L tot - L core ), where L tot is the total (single-dish) 
luminosity. We list the logarithm of this quantity (in the 
source frame) in column 10 of Table 2. 

A beneficial aspect of the PR sample is that high-quality 
arcsecond-scalc 1.4 GHz images have been obtained for 
nearly all of the sources as part of the Caltech-Jodrcll sur- 
veys (Polatidis et al. 1995) or other studies. In columns 
11 and 12 of Table 2 we list the 1.4 GHz luminosities of 
the arcsecond-scale core and extended structure, respec- 
tively, as measured from VLA or MERLIN images. The 
majority of these data are from Murphy et al. (1993), who 
obtained high dynamic range images of a large sample of 
core-dominated AGNs with the VLA. Here we assume typ- 
ical spectral indices of a = for the core components, and 
a = —0.8 for the extended emission. For those objects 
with no detectable extended emission, we assumed an up- 
per limit of three times the rms noise level in the map. This 
limit is valid provided the undetected emission is spread 
over an area larger than the restoring beam. 

In column 13 we tabulate the kiloparsec-scale core dom- 
inance ratios (Rk pc ) for our objects, which is simply the 
ratio of core to extended luminosity at 1.4 GHz, measured 
in the rest frame of the source. In column 15 we list the 
monochromatic luminosities of our sources in the 1 keV X- 
ray band. These were derived from fluxes in the literature, 
assuming a xray — —0.8. 

We have also calculated optical luminosities (not listed) 
based on optical V band magnitudes from NED 6 . We made 
no correction for galactic absorption, and adopted a stan- 
dard photometric transformation (Wamsteker 1981) and 
an optical spectral index a — —0.5. 

3.5. Maximum apparent jet velocity 

The apparent velocities of distinct components in AGN 
jets, combined with estimates of the viewing angle, can 
provide useful information on the bulk Lorentz factor of 
the flow. Many superluminal AGNs (e.g., 3C 345; Ros, 

5 http://www.astro. lsa.umich.edu/obs/radiotcl/umrao. html 
6 http://ned www. ipac.caltech.edu/ 



Zensus, & Lobanov 2000) exhibit components moving at 
different velocities, which suggests that we are likely seeing 
a pattern speed rather than the true speed of the under- 
lying plasma. Nevertheless, the maximum observed veloc- 
ity can still provide a useful constraint on the amount of 
relativistic beaming in a jet. Apparent velocity data are 
available for 23 sources in the VSOP-PR and FS-PR sam- 
ples. In column 2 of Table 3 we list the maximum velocity 
for these sources in units of c. 

3.6. Jet bending 

There is a great deal of interest in mapping the apparent 
trajectories of AGN jets at high resolution, as small varia- 
tions in viewing angle can dramatically alter the apparent 
brightness distribution along a relativistic jet (Gomez et 
al. 1994). Many authors (e.g., Pearson & Readhead 1988; 
Tingay et al. 1998) have used APA, the difference be- 
tween the jet position angles measured on milliarcsccond 
and arcsecond scales, as a simple measure of jet bending. 
For those sources that have either straight jets or a single 
bright component on arcsecond scales, the measurement 
of APA is straightforward. However, as pointed out by 
Piner (1998), the wide variety of complex curvatures and 
morphologies found in AGN jets makes it difficult to mea- 
sure APA in a uniform manner for large samples. For 
some curved jets such as 1633+382, there have been up 
to five different values of APA published in the literature. 
We have attempted to limit possible biases due to redshift 
in our analysis by defining the position angle on kiloparsec 
scales (PAkpc) to be the position angle of the arcsecond- 
scale jet at a projected distance of 10 hr 1 kpc from the 
core. In cases where there was no identifiable jet, we used 
the position angle of the first radio component located far- 
ther than this projected distance. We believe this is the 
most self-consistent method, given the large range of red- 
shift range (and in turn, image surface brightness) of our 
sample. 

We define the position angle on parsec scales (PA pc ), 
as the position angle of the innermost jet component with 
respect to the core. Since the region near the core in most 
of our sources is still optically thick at 5 GHz, we are cur- 
rently obtaining 43 GHz images of the FS-PR sample with 
the VLBA to obtain better estimates of PA pc (Lister et 
al., in preparation). In column 2 of Table 3 we list the 
innermost jet position angle at 43 GHz for those sources 
which we already have data, otherwise we give a position 
angle based on our 5 GHz VSOP image or other published 
VLBI data. We list the APA values for our sources in 
column 8 of Table 3. 

3.7. Polarization properties 

Extensive data on the linear polarization properties of 
the PR sample in the optical and radio regimes have been 
tabulated by Impey et al. (1991) and Aller et al. (1992), 
respectively. In Table 1 we list the measured properties 
from these papers that we include in our correlation anal- 
ysis. Our optical polarization data are taken from Impey 
et al. (1991) with the exception of 1954+513, whose data 
are from Wills et al. (1992). We introduce the derived 
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quantity |%15 — PA pc \, which represents the angular offset 
between the inner jet position angle and the electric po- 
larization vector angle (x) at 15 GHz. The X15 values are 
time-averaged single-dish quantities tabulated by Aller et 
al. (1992). Due to the 180° ambiguity in the x measure- 
ments, \xi5 — PA pc \ cannot exceed 90°. 

4. CORRELATION ANALYSIS 

4.1. Description of statistical tests 

The large variety of data that we have gathered on 
the FS-PR sample are well suited for a multi-dimensional 
correlation analysis. This method checks for correlations 
among all pairwise combinations of variables, and is there- 
fore unbiased with respect to any a priori expectations. A 
major caveat, however, is the possibility of detecting spu- 
rious correlations. For example, when testing for correla- 
tions among n variables, if one chooses a significance level 
p > P to call a significant correlation, you should expect 
to find (1 — P)n(n — l)/2 spurious correlations. Unfortu- 
nately, it is usually impossible to determine which corre- 
lations are spurious, and which arc genuine. If the dataset 
is large enough, it is sometimes feasible to run correlation 
tests on different subsets of the data. Alternatively, the 
tests can be run on an entirely new dataset. Neither of 
these options is available to us for the FS-PR sample, so 
we have reduced the chance of false detections by adopting 
a relatively high confidence level of 98%. For every corre- 
lation above this significance, we also checked the scatter 
plot to determine whether the significance level had been 
artificially increased by the presence of an outlying data 
point. If the significance dropped below (98%) after the 
removal of the outlier, we rejected the correlation. 

To determine the correlation significance among pairs of 
variables, we employed a Kendall's tau test. This test does 
not require that the intrinsic relation between variables be 
linear (i.e., it is non-parametric), and considers only the 
relative ranks of individual data pairs. Since many of our 
data points are censored (i.e., are limits only), we cm- 
ployed a special version of Kendall's tau that incorporates 
the technique of survival analysis. We used the algorithms 
provided in the ASURV package (Lavalley et al. 1992) for 
this purpose. 

Many of the measured properties of flux-limited AGN 
samples have a large dependence of redshift, due to the 
steepness of the parent luminosity function. Intrinsically 
faint radio sources in these samples typically have low red- 
shifts due to the flux cutoff. As as result, certain variables 
with a mutual dependence on redshift, such as luminosi- 
ties at different wavelengths, can appear to be very well- 
correlated, even though no intrinsic correlation may exist 
(see, e.g., Padovani 1992). In order to eliminate these 
rcdshift-induced correlations from our analysis, we calcu- 
lated partial correlation coefficients of the form 

7~ab TazTfrz . , 

T ab,z — , ■ , , \£) 

where the r's represent the various Kendall's tau corre- 
lation coefficients between variables a,b, and redshift (z). 
Extra steps must be taken when calculating partial corre- 
lations involving censored data. We used the algorithms 
of Akritas & Siebert (1996) that have been developed for 
this purpose. 



In addition to our Kendall's tau tests, we also performed 
statistical tests to look for differences in source proper- 
ties among different classes, e.g., BL Lacs versus quasars 
and EGRET versus non-EGRET sources. For those vari- 
ables that involved limits, we used Gehan's generalized 
Wilcoxon two-sample test (Gchan 1965) from the ASURV 
package. For the remaining variables we used a standard 
Kolmogorov-Smirnov (K-S) test. In both cases we adopted 
a 98% confidence level in determining whether the differ- 
ences were significant. 

4.2. Results 

Our Kendall's tau tests have identified 30 out of 171 
possible variable pairs that are correlated above the 98% 
confidence level (we discuss additional correlations with 
redshift separately in §4.2.1). From our previous discus- 
sion in §4.1, we can expect approximately three of these 
to be due to chance. Of the original 30 correlations, we 
reject six as likely being due to a mutual dependence on 
redshift, and another three as a result of one or two out- 
lier points that artificially boost the confidence level. This 
leaves 21 bone fide correlations for the sample, 20 of which 
are either expected or have been previously found in other 
studies (see Table 4). 

We have resorted to several different methods in order 
to summarize the complex inter-relations among this large 
collection of source properties. Since a single correlation 
matrix is too large to present here, we list in Table 4 only 
the 30 initial pairs that are significant at the 98% level or 
higher. We give a reference for each correlation that has 
been previously detected for radio-loud AGNs by other au- 
thors. We note that in some cases this may not refer to 
the very first detection of a particular correlation. 

In Table 5 we summarize the results of our two-sample 
tests. The values represent the probabilities (in per cent) 
that the samples were drawn from the same parent pop- 
ulation. For clarity, we have omitted all values from this 
table that exceed 2%. 

In Figure 1 we present a graphical representation of our 
test results. The solid lines represent correlations found 
from Kendall's tau tests, while the broken lines denote 
those from the two-sample tests. Arrows on the lines 
denote anti-correlations, and the lines in bold represent 
new correlations detected in this study. For example, the 
dashed line connecting IDV and m opt indicates that IDV 
type I and II sources have higher optical polarization lev- 
els than type sources. There are indications of several 
groupings of variables in this diagram, which we will dis- 
cuss in more detail in §5.1. 

4.2.1. Correlations with redshift 

All of the luminosity variables in our dataset, with the 
exception of extended luminosity (L ext ), are correlated 
with redshift at better than 99.999% confidence. This 
is expected, given the flux-limited nature of the FS-PR 
sample. None of the other (non-luminosity) variables are 
significantly correlated with redshift. We note that the BL 
Lacertae objects have significantly lower redshifts than the 
quasars in the sample (Table 5). 

4.2.2. Correlations involving luminosity variables 
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Fig. 1. — Diagram summarizing interrelations among source properties for the FS-PR sample (see Table 1 for symbol definitions). Solid 
lines indicate correlations found from Kendall's tau tests. Dashed lines indicate significant population differences according to two-sample 
tests. Arrows denote anti-correlations, and lines in bold represent new correlations detected in this study. 
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The various luminosities of our sample objects are gen- 
erally well-correlated with each other, but not with other 
source properties, as is indicated by their location in the 
lower right-hand corner of Figure 1. One major exception 
is IDV activity. Our two-sample tests indicate that sources 
that have not shown IDV tend to be more luminous at 320 
MHz and 5 GHz. These trends may be due to beaming 
and selection effects in the sample (see §5.3). 

Among luminosity-luminosity pairs, correlations involv- 
ing the total 5 GHz luminosity (L5) are the most prevalent. 
With the exception of i 5 versus the VLBI core luminosity 
{L c , P c), ah of the luminosity- luminosity correlations in Ta- 
ble 4 have been previously detected in other studies. The 
strong correlation between L5 and L c pc is merely a re- 
flection of the highly core-dominated nature of the FS-PR 
sample. 



4.2.3. Trends with polarization 

The level of polarization at optical wavelengths has his- 
torically been an effective method of identifying highly 
beamed blazars whose jets are seen nearly end-on, since 
very few optically-selected quasars have m opt > 3% 
(Stockman et al. 1984). In Table 4 we confirm the pre- 
viously known correlations with parsec-scale core domi- 
nance, 5 GHz radio variability amplitude, and emission 
line equivalent width that have helped establish m opt as a 
beaming indicator (see, e.g., Impey et al. 1991). 

Xu et al. (1994) found that in the combined 
PR/Caltech-Jodrell samples, sources with low optical po- 
larization levels (< 3%) tend to have better jet alignments 
between parsec and kiloparsec scales than highly polar- 
ized sources. We find some indication of differences of 
these groups in the m opt - APA plane (Fig. 2). There is 
a distinct grouping of sources in the lower left-hand por- 
tion of the plane which have both low optical polarization 
and small jet misalignments. The remainder of the sample 
shows a roughly linear trend of decreasing optical polar- 
ization with increasing jet misalignment. 

It is unclear whether the differences between the high 
and low-optically polarized objects in Fig. 2 can be entirely 
attributed to beaming and projection effects, which would 
tend to amplify any small intrinsic jet bends in sources 
seen nearly end-on. The current consensus is that blazars 
appear highly polarized since their optical synchrotron ra- 
diation is significantly boosted with respect to the underly- 
ing (thermal) continuum (the same process is responsible 
for lowering the observed emission line widths of highly 
beamed objects, e.g., Wills et al. 1992). If the optical syn- 
chrotron component is not highly polarized to begin with, 
however, the optical polarization will remain low, regard- 
less of how beamed the source is. For example, Lister & 
Smith (2000) found evidence that many of the differences 
between HPQs and LPRQs are not due to differences in 
Dopplcr factor or jet orientation, but rather intrinsic dif- 
ferences in their jet magnetic field structures. It is pos- 
sible that blazars go through temporary quiescent phases 
of low intrinsic polarization and flux variability, in which 
they appear as LPRQs (e.g., Fugmann 1988). Long-term 
optical polarization monitoring studies of samples such as 
the FS-PR are needed to establish the properties of this 
duty cycle in blazars. 



4.2.4. Properties of IDV sources 

In an accompanying paper (Tingay et al. 2001), we dis- 
cuss how the PR sources that have displayed type I and II 
IDV have substantially higher VSOP core brightness tem- 
peratures than those sources that have not displayed any 
IDV. Although it appears that the presence of a bright, 
compact core on VLBI scales is a necessary criterion for 
IDV (e.g., Heeschen et al. 1987; Quirrenbach et al. 1992), 
the phenomenon only occurs in <; 2/3 of all flat-spectrum 
radio sources (Quirrenbach et al. 1992; Kraus et al. 1999). 
Based on the IDV statistics of a large sample of AGNs, 
Quirrenbach et al. (1992) concluded that IDV was more 
prevalent in sources with 80% or more of their total 5 GHz 
flux density contained in an unresolved VLBI core. We 
detect a similar trend in the FS-PR, with the R pc distri- 
bution of the type (non-IDV) sources differing from that 
of the type II sources at the 99.918% confidence level. The 
type sources also have steeper radio spectral indices, in 
agreement with the findings of Quirrenbach et al. (1992). 

We have also detected two additional trends that con- 
firm the association of IDV with the blazar class. The IDV 
type I and II sources have higher optical polarization lev- 
els and narrower emission line equivalent widths than the 
type objects. This suggests that the weak-lined BL Lac 
objects may be the best candidates for IDV; indeed, an 
optical study by Heidt & Wagner (1996) found IDV to be 
present in 28 of 34 BL Lacs selected from the 1 Jy catalog 
(Kiihr et al. 1981). 

The correlations with optical polarization and equiva- 
lent width suggest that the occurrence of IDV depends 
more on the Doppler factor of the source rather than any 
intrinsic property such as core size. As we will show in 
§5.3, the highest Doppler factor sources in a flux-limited 
sample such as the FS-PR are predicted to have low intrin- 
sic (unbeamed) luminosities, according to the relativistic 
beaming model. Since the 320 MHz luminosity is gener- 
ally accepted to be a good indicator of total jet power and 
intrinsic luminosity (Rawlings & Saunders 1991; Kaiser & 
Alexander 1997), this would imply that the type I and II 
IDV sources in the FS-PR should have lower L 32 o values, 
which is in fact the case (Fig. 3). 

It is significant that we found no trends between IDV 
and beaming indicators on kiloparsec scales (i.e., APA, 
Rkpc, arid morphological classification). This further im- 
plies that it is the Doppler factor of the innermost region 
of the jet that determines whether a source will display 
IDV. 

4.2.5. Dependence of source properties on extended radio 

structure 

Our correlation analysis has revealed that the arcsecond- 
scale radio morphology of core-dominated radio sources 
is closely connected with several other important source 
properties. In §3.1 we classified our sources as either 
compact, one-sided, two-sided, or halo based on their ex- 
tended arcsecond-scale radio structure. Murphy (1988) 
analyzed the kiloparsec-scale properties of a large sample 
of core-dominated AGNs, and found that the two-sided 
sources tend to have lower core-dominance ratios {Rkpc) 
and stronger extended luminosities. We confirm these 
findings for the FS-PR, and additionally find that the two- 
sided sources have smaller parsec-scale core dominance, 
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larger emission line equivalent widths, and smaller jet mis- 
alignments between parsec- and kiloparsec-scales (Fig. 4) . 
Impey et al. (1991) found that the latter three proper- 
ties are all highly influenced by source orientation, which 
implies that the jet axes of the two-sided sources are ori- 
ented farther away from the line of sight than the other 
sources in the FS-PR. It is conceivable that the two-sided 
sources might have intrinsically smaller core dominance 
and stronger extended power (irrespective of beaming), 
but we consider it unlikely that intrinsic differences could 
also conspire to mimic the trends with equivalent width 
and jet misalignment. Furthermore, there is strong evi- 
dence from unbeamed samples such as the 3C survey that 
the vast majority of powerful radio loud AGNs are intrin- 
sically two-sided (Hintzen, Ulvestad & Owen 1983). A 
simple explanation for the trends we see with sidedness is 
that for sources with jets seen nearly end-on, the emission 
from the hotspot on the far side of the one-sided sources is 
either beamed away from us, or is swamped by the highly 
beamed core in images of limited dynamic range. 

4.2.6. Jet misalignment and bending 

The difference in inner jet position angle between parsec 
and kiloparsec scales (APA) has been used in many stud- 
ies as a general indicator of jet orientation. Since AGN 
jets are known to be intrinsically bent, those jets that 
are closer to the line of sight should have more exagger- 
ated bends due to projection effects, and usually (but not 
necessarily) larger values of APA. Pearson & Readhcad 
(1988) found unexpectedly that the distribution of APA 
for the PR survey was bi-modal, with a clear division at 
approximately 45 degrees. Xu et al. (1994) found a sim- 
ilar result for the combined PR/Caltech-Jodrell samples. 
The reasons for this bi-modality are still unclear. Conway 
& Murphy (1993) have suggested that it could be due to 
selection effects associated with relativistic beaming and 
jet trajectories in the form of damped helices. 

In the top panel of Figure 2 we show the distribu- 
tion of APA for the FS-PR sample, based on our im- 
proved parsec- and kiloparsec-scale jet position angle mea- 
surements described in §3.6. The bi-modality noticed by 
Pearson & Rcadhead (1988) is still present, with a divi- 
sion between aligned and misaligned sources occurring at 
roughly 60°. The vast majority of the aligned population 
(APA < 60°) have two-sided radio structure on arcsecond- 
scales (Fig. 4). This dependence of jet misalignment on 
lobe morphology strongly suggests that the highly mis- 
aligned sources have inner jets that are oriented closer to 
the line of sight. The observed trend of APA with emis- 
sion line equivalent width (also noticed by Appl et al. 1996 
in a larger AGN sample) also supports this interpretation. 

4.2.7. Parsec-scale core dominance 

Numerous studies (e.g., Impey et al. 1991; Wills et al. 
1992) have found that parsec-scale core dominance (R pc ) 
is likely a good beaming indicator, given that it is well- 
correlated with several other source properties that are af- 
fected by beaming. These include a strong anti-correlation 
with emission line equivalent width (Impey et al. 1991), 
and a positive correlation with optical polarization (Im- 
pey et al. 1991), both of which we confirm here (Table 4). 
In Figure 5 we show the distribution of parsec-scale core 



dominance for the various optical classes in the FS-PR 
sample. Our K-S tests reveal no statistically significant 
differences in these distributions. 

4.2.8. Kiloparsec-scale core dominance 

In Figure 6 we show the distribution of kiloparsec-scale 
core dominance for the FS-PR. The measured values and 
limits span a remarkable range of at least four orders of 
magnitude, which reflects the wide variety of extended 
structure present in the sample (the parsec-scale core dom- 
inance spans only ^2.5 orders of magnitude in compari- 
son). We confirm the finding of Morganti et al. (1995) that 
more highly core-dominated sources tend to have weaker 
extended luminosities. 

The core dominance on kiloparsec scales is not corre- 
lated with any parsec-scale quantity for the FS-PR, such 
as emission line equivalent width, or optical polarization, 
as is the case for the parsec-scale core dominance. These 
correlations have been found, however, for non-core se- 
lected samples, e.g., Browne & Murphy (1987). It appears 
that VLA core dominance is not a useful indicator of in- 
ner jet orientation for samples selected on the basis of core 
flux density. 

4.2.9. Apparent jet velocity 

The maximum apparent jet velocity of a source is not 
necessarily a good indicator of relativistic beaming and jet 
orientation in core-selected samples, due to the non-linear 
dependence of apparent velocity (P apP ) on viewing angle. 
Since the apparent velocity varies as (3 cos 9/(1 — /3cos9), 
where (3 is the true jet velocity in units of c and 9 
is the viewing angle, highly beamed jets aligned within 
1/r = a/1 — ft 2 radians to the line of sight can have rather 
low apparent velocities. We did not detect any correlations 
between apparent velocity and any other source property 
for the 20 FS-PR sources for which proper motion data 
are available. It appears that the full PR sample, which 
contains a wider range of jet viewing angles, may be better 
suited for exploring trends with apparent velocity. 

4.2.10. Gamma-ray emission 

The FS-PR sample contains seven AGNs that were de- 
tected in high-energy gamma-rays by the EGRET obser- 
vatory. Our two-sample tests showed no significant differ- 
ences in the properties of the EGRET versus non-EGRET 
objects. It is possible, however, that this null result may be 
due to the sample sizes (7 versus 25 objects), since the sta- 
tistical reliability of Kolmogorov-Smirnov tests is greatly 
decreased for such small samples (Press et al. 1992). 

5. DISCUSSION 

5.1. Cluster analysis 

Our network diagram (Fig. 1) graphically displays the 
relationships between the various properties of the FS-PR 
sources. In order to investigate possible groupings in this 
diagram, we have applied a cluster analysis to our data 
(e.g., Nair 1997). This method searches for natural group- 
ings among variables that are most correlated. The most 
common ( "single- link" ) method involves taking a symmet- 
ric matrix of correlation coefficients and finding the two 
variables a and b that are the best correlated. These two 
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variables are replaced by a new "cluster" variable c, for 
which r ci = max[r oi , r^]. The process is repeated until 
there is only one cluster variable remaining. The results 
are usually presented in the form of a tree diagram, in 
which the clusters are represented by branches. 

A few modifications to the standard cluster method were 
needed to apply it to our data. Since our computed cor- 
relation coefficients for censored variables have a differ- 
ent form of variance than the uncensored ones (Akritas 
& Siebert 1996), they cannot be directly inter-compared. 
Also, the fact that many of our variables have a mutual 
dependence on redshift will artificially skew the shape of 
the cluster tree diagram. We have therefore taken the 
following steps in order to arrive at a more accurate rep- 
resentation of the true cluster tree. First, we have used 
the statistical significance levels in place of the correlation 
coefficients to compare the variables. For the two nom- 
inal variables (IDV and extended morphology class), we 
have used the two-sample test significance levels. Second, 
we have used the significance levels with redshift partialcd 
out for all of the luminosity variables. Finally, we have 
replaced the significance level for the rejected correlations 
listed in Table 4 with a value that excludes the outlier 
point. 

We show the resulting tree diagram in Figure 7. The 
vertical connections indicate a close connection between 
variables and groups of variables. The horizontal axis is 
roughly proportional to correlation strength, with the ver- 
tical bars located furthest to the right representing the 
strongest correlations. The luminosity variables form a 
distinct grouping since they are poorly correlated with 
the remaining variables. The apparent jet velocity and 
fractional polarization levels at 5 & 15 GHz also form a 
distinct group that is largely independent of other source 
properties. 

Among the remaining variables, those associated with 
the parsec-scale core form a distinct grouping (labeled 
"parsec scale" in Fig. 7). The emission line equivalent 
width, kiloparsec-scale morphology, and jet misalignment 



angle form a separate grouping. The exclusion of the 
equivalent width from the parsec-scale grouping is some- 
what unexpected, given that the broad line region in 
quasars is thought to be less than a parsec in extent (Capri- 
otti, Foltz & Peterson 1982). This may imply that the 
kiloparsec-scale properties of jets such as bending are in- 
fluenced by the intrinsic line strength, and in turn the 
density of the broad line gas. 

5.2. Emission line strength and the BL Lacertae class 

Early VLBI studies of AGN samples with complete op- 
tical spectra (e.g. Impey et al. 1991) showed that emis- 
sion line strength was correlated with core dominance on 
parsec-scales, and led to a model in which the optical con- 
tinuum in highly beamed sources is dominated by syn- 
chrotron emission from the jet (Wills et al. 1992). Since 
the emission from the broad line region is unbeamed, the 
ratio of line to continuum flux (as well as emission line 
equivalent width) should steadily decrease with viewing 
angle to the source. Also, the overall optical polarization 
should increase, as the polarized synchrotron component 
begins to dominate. All of these predictions are borne out 
in the FS-PR, as evidenced by the mutual correlations be- 
tween equivalent width, APA and extended morphology. 

These correlations have potentially important implica- 
tions for the BL Lacertae class of AGN, which by definition 
have low emission line equivalent widths. Several studies 
have shown that these objects have intrinsically different 
properties than their broad- lined counterparts, the high 
optically-polarized quasars (HPQs). These include differ- 
ences in their parsec-scale magnetic fields (Cawthorne et 
al. 1993), X-ray spectral shapes (Worrall & Wilkes 1990), 
and extended radio power (Padovani 1992). It has been 
suggested that the parent population of the BL Lacs are 
FR-I type radio galaxies, while the radio quasars belong to 
the more powerful FR-II population (see Urry & Padovani 
1995 and references therein). 

The nine BL Lacs in the FS-PR rank in the top 30% of 
the sample in terms of optical polarization, APA values, 
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5 GHz variability amplitude, and parsec-scale core domi- 
nance. This implies that they are indeed highly beamed 
objects. There remains some debate over whether BL Lacs 
are in fact more beamed (i.e., have higher Doppler factors) 
than high-optically polarized quasars. The FS-PR sample 
is not ideally suited to fully investigate this issue, since 
it contains only nine BL Lacs and nine confirmed HPQs. 
Nevertheless, we have performed several K-S tests on our 
data, and find that although the BL Lacs do tend to have 
lower than average redshifts, there are otherwise no sig- 
nificant differences in the APA, optical polarization, vari- 
ability amplitude, and core dominance distributions of the 
two classes. Therefore, despite strong indications that the 
observed equivalent width in core-dominated radio sources 
is highly dependent on beaming, there is insufficient evi- 
dence in the FS-PR to support a claim that BL Lacs are 
more highly beamed than HPQs. This hypothesis should 
not be ruled out, however, until a similar analysis can be 
conducted on a much larger sample such as the Caltech- 
Jodrcll Flat-spectrum survey (CJ-F; Taylor et al. 1996). 

5.3. Relativistic beaming simulations 

In the previous discussion we have made reference to 
general predictions of the relativistic beaming model, and 
how it can account for many of the correlations seen in 
the FS-PR sample. In this section we describe in more 
detail the predicted properties of sources in core-selected, 
flux-limited samples, based on Monte Carlo beaming sim- 
ulations. 

Simulations of relativistic beaming by Vermeulen & Co- 
hen (1994) and others have shown that core-selected sam- 
ples should be dominated by objects with jets pointing 
nearly directly at us. This would suggest that the FS-PR 
should contain a narrow range of highly beamed, high- 
Doppler factor sources. It is essential to consider, however, 
the additional effects of the luminosity function (LF) and 
rcdshift distribution of the FS-PR parent population. It is 
very likely that some low-Doppler factor sources that hap- 
pen to be intrinsically very luminous or relatively nearby 
may have also been included in the FS-PR sample. This 
would explain the presence of the two radio galaxies (3C 
84 and 2021+614), which are generally considered to be a 
weakly beamed class of AGN (Urry & Padovani 1995). It 
would also account for the large range of kiloparsec-scale 
core dominance found within the sample. The effects of 
adding a LF and redshift distribution to beaming models 
have been investigated by Lister & Marscher (1997). They 
found that samples such as the FS-PR should in fact con- 
tain a wide range of Doppler factors. In this section we 
use the model of Lister & Marscher (1997) to show how 
many of the correlations in the FS-PR are likely a result 
of a mutual dependence of source properties on Doppler 
factor. 

Lister & Marscher (1997) used their Monte Carlo model 
to successfully fit various properties of the Caltech-Jodrcll 
Flat-spectrum sample (CJF; Taylor et al. 1996), of which 
the FS-PR is a sub-set. Their best-fit parent population 
contains a power law jet Lorcntz factor distribution of the 
form N(T) oc r^ 125 with 1.001252 < T < 30. The jets 
are randomly oriented in space, and distributed with a 
constant co-moving space density out to z = 4. The jet 
luminosity function mirrors that of the FR-II radio galax- 



ies, which was derived by Urry & Padovani (1995), and 
incorporates pure exponential luminosity evolution. The 
jets are assumed to be two-sided, and have luminosities 
that are Doppler boosted by a factor of S 2 . In a typical 
model run, parent objects are continuously created until 
293 sources with flux density exceeding 0.35 Jy (the CJF 
sample size and cutoff) are obtained. Further details of 
the model can be found in Lister & Marscher (1997). 

In Figures 8, 9, and 10 we plot the results of a single 
Monte Carlo simulation with the above parameters that 
illustrates several of the strong biases that are present in 
flux-limited samples. The small circles represent the 293 
simulated CJF sources, and larger circles indicated sources 
brighter than 1.3 Jy (the FS-PR cutoff). Figure 8 plots 
the intrinsic luminosity of the jet (removing the Doppler 
boosting effect) against Doppler factor. The upper enve- 
lope to the distribution indicates that extremely powerful 
sources that are also highly beamed are exceedingly rare, 
and that weakly beamed sources generally need to be more 
intrinsically luminous to be included in flux-limited sam- 
ples. This trend provides a natural explanation for the ob- 
servation that IDV type sources generally have stronger 
320 MHz and 5 GHz luminosities than the IDV type I and 
II sources in the FS-PR, since the non-IDVs are expected 
to be less relativistically boosted (see §4.2.4). 

In Figure 9 we plot two quantities against the jet Lorentz 
factor: apparent velocity (upper panel), and viewing an- 
gle (lower panel). The solid lines on these panels indicate 
the positions of sources seen at the viewing angle that 
maximizes the apparent velocity (9 = arcsin [1/F]). As 
pointed out by Lister & Marscher (1997), there is a sub- 
stantial spread around this angle (lower panel), with most 
sources seen at significantly smaller angles to the line of 
sight. Nevertheless, the fairly good correlation in the up- 
per panel suggests that apparent velocity can be used as a 
statistical indicator of the Lorentz factor. We note, how- 
ever, that our simulations assume the observed pattern 
speed reflects the true speed of the jet, which may not be 
the case for all sources (Vermeulen & Cohen 1994). 

The apparent velocity of a jet by itself is not a good 
indicator of viewing angle (lower panel of Fig. 10), due 
to the preponderance of jets seen within the 1/r view- 
ing cone. At best it can be used to set an upper limit 
on the viewing angle, through the well known formula 
9 < arccos [{Pl PP - 1)/ {Pa PP + !)]• The lack of a good cor- 
relation between apparent velocity and viewing angle may 
explain why the former quantity is not correlated with 
other source properties in the FS-PR. Alternatively, the 
measured apparent speeds could simply be pattern speeds 
that do not reflect the true jet Lorentz factor. If this is the 
case, the best statistical approach would be to use long- 
term VLBI monitoring to determine a maximum observed 
velocity for each source in a large sample. 

The upper panel of Fig. 10 shows that we can ex- 
pect an extremely good correlation between viewing an- 
gle and Doppler factor in a small flux-limited sample such 
as the FS-PR, provided the parent population has an up- 
per Lorentz factor limit. This trend is likely the source 
of correlations between orientation-dependent properties 
(e.g., jet misalignment angle and arcsecond-scale mor- 
phology) and Doppler boosting-dependent properties (e.g., 
core dominance and emission line equivalent width). 
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5.4. The identification of highly beamed radio sources 

The Doppler factors of extragalactic jets are notoriously 
difficult to measure accurately, due to a lack of good indi- 
cators of true jet velocity and orientation. Although some 
highly beamed sources can be recognized on the basis of 
blazar properties (i.e., high variability and optical polar- 
ization), other beamed sources may have intrinsically low 
variability and/or polarization. Our analysis of the FS-PR 
sample suggests a set of criteria that can be used to identify 
the AGNs with the highest Doppler factors in core-selected 
samples. In general, the most highly beamed sources will 
tend to have a) flat radio spectra, b) a type I or II IDV clas- 
sification; c) one-sided, halo-like, or no extended structure 
on kiloparsec scales; d) a relatively large (;> 60°) misalign- 
ment between their jet directions on parsec and kiloparsec 
scales; e) optical polarization that has exceeded ~ 3% on 
at least one occasion; f) relatively high core dominance on 
parsec (R pc ^ 0.3) and kiloparsec (Rk pc ^ 10) scales; g) 
an equivalent width of the largest permitted emission line 
smaller than ~ 15 A (source frame); and h) a 5 GHz source 
frame core brightness temperature exceeding 10 11 K. 

After taking into account the missing data on some of 
our sources, there are three HPQs (0133+476, 0804+499, 
and 1739+522) and six BL Lacs (0454+844, 0954+658, 
MK 501, 1749+701, 1803+784, and BL Lacertae) in the 
FS-PR that are consistent with the above criteria. We sug- 
gest that these sources have significantly higher Doppler 
factors than the other sources in the FS-PR, and should 
be considered as prime targets for high-frequency space- 
VLBI and IDV monitoring observations. We predict that 
two of these objects (0454+844 and 1739+522) will even- 
tually be found to display type I or II IDV, based on their 
other properties. 

6. CONCLUSIONS AND SUGGESTIONS FOR FUTURE 
WORK 

We have obtained data using the facilities of the VSOP 
mission and from the literature on a flat-spectrum, core- 
selected subsample of the Pearson-Readhead AGN survey 
(the FS-PR sample) in order to investigate the effects of 
relativistic beaming in compact radio sources. The main 
results of our multi-dimensional correlation analysis on 
these data are as follows: 

(i) The properties of sources that have displayed type I 
or II IDV are significantly different than those that have 
shown no IDV activity. They have flatter spectral indices, 
higher core dominance on parsec scales, higher optical po- 
larizations and smaller emission line widths. They also 
have lower 320 MHz and 5 GHz luminosities and weaker 
kiloparsec-scale core components. These findings are all 
consistent with Monte Carlo beaming simulations which 
predict that higher Doppler factor sources in flux-limited 
generally have weaker intrinsic luminosities. 

(ii) The observed luminosity properties of the FS-PR 
are generally well-correlated with each other, but not with 
other properties of the sample. The only exception is the 
trend with IDV described above. 

(iii) The high-optically polarized sources in the FS-PR 
show a trend of increased parsec-to-kiloparsec jet misalign- 
ment (APA) with decreasing optical polarization. How- 
ever, the lowest optically polarized sources (m op t < 3%) 
generally have small jet misalignments, which may suggest 



that they have intrinsically different jet properties. 

(iv) We find that sources with two-sided radio morphol- 
ogy on kiloparsec-scales have significantly different prop- 
erties than other sources in the FS-PR. They tend to be 
less core dominated on parsec- and kiloparsec scales, and 
have stronger extended luminosities and larger emission 
line widths. Their jets are also less misaligned between 
parsec- and kiloparsec scales, which suggests that they are 
oriented at angles farther away from the line of sight. 

(v) We have used high-resolution VLBI images to con- 
firm the bi-modality of the APA distribution for the PR 
sample (Pearson & Readhead 1988), with a division oc- 
curring APA ~ 60°. The aligned sources tend to have 
broader emission line equivalent widths than the mis- 
aligned population. 

(vi) We did not find any correlations between maximum 
apparent jet velocity and other source properties. This 
may be due to the predicted weak dependence of (3 app on 
viewing angle for core-selected samples, or the fact that 
the measured speeds may not reflect the true bulk speeds 
of the emitting material in the jet. 

(vii) The nine BL Lacertae objects in the FS-PR rank 
among the top 30% of the sample in terms of their opti- 
cal polarization, APA values, 5 GHz variability ampli- 
tude, and parsec-scale core dominance, which suggests 
that they are highly beamed objects. We found no sig- 
nificant differences in their properties as compared to the 
high-polarization quasars, other than the fact that they 
have smaller redshifts. 

(viii) Our correlation analysis and Monte Carlo beaming 
simulations have shown that the properties of the FS-PR 
sample are remarkably consistent with the predictions of 
the relativistic beaming model. The majority of the ob- 
served correlations are likely due to a mutual dependence 
of source properties on the jet Doppler factor and view- 
ing angle, and do not necessarily reflect intrinsic relations 
among jet properties. 

Our work has demonstrated the merits of obtaining a 
wide variety of data on a relatively small sample, and sug- 
gests several further observations that can improve our un- 
derstanding of the beaming phenomenon in AGNs. These 
include optical variability monitoring of the PR survey to 
further investigate the connection between optical and ra- 
dio properties, and more uniform monitoring of the sample 
for IDV activity in the radio. It will also be important to 
measure apparent velocities for the remaining sources in 
the sample, and to obtain high-dynamic range VLA and 
MERLIN maps down to a uniform luminosity level. The 
latter can be used to trace out bending in the jets as they 
propagate from parsec- to kiloparsec scales, and examine 
how this phenomenon is influenced by other source prop- 
erties. With regular monitoring the optical polarization of 
the sample at regular intervals, it should be possible to es- 
tablish a duty cycle of optical polarization in blazars, and 
to better understand the inherent differences in high- and 
low-optically polarized AGNs. 
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Table 1 

Measured and Derived Quantities for the Pearson-Readhead AGN Sample 



Symbol 



Property 



Units 



Reference 



General Properties 

EGRET gamma-ray detection/non-detection 
Optical classification (RG, LPRQ, HPQ, BL Lac) 
Kiloparsec-scale extended morphology 
z Redshift 

EW Equivalent width of widest permitted line (source frame) 

a 5 _i 5 Time-averaged spectral index between 4.8 and 14.5 GHz 

Variability Properties 
V 5 Variability amplitude at 5 GHz 

V15 Variability amplitude at 14.5 GHz 

IDV classification 

Luminosity Properties 



A 



Hartman et al. 1999 

Various 

Various 

Lawrence et al. 1996 
Lawrence et al. 1996 
Allcr et al. 1992 



Aller et al. 1992 
Aller et al. 1992 
Various 



£320 


Total luminosity at 320 MHz 


WHz" 


-1 


Impey et al. 1991 


L 5 


Total luminosity at 5 GHz 


WHz" 


-1 


UMRAO database a 


L pt 


Optical V band luminosity 


WHz" 


-1 


NED 7 


Lxray 


X-ray luminosity at 1 keV 


WHz" 


-1 


Various 


Lc.pc 


Parsec-scale core luminosity at 5 GHz 


WHz" 


-1 


This paper 


Lc.kpc 


Kiloparsec-scale core luminosity at 1.4 GHz 


WHz" 


-1 


Various 


Lext 


Total arcsecond-scale extended luminosity at 1.4 GHz 


WHz" 


-1 


Various 



Jet Properties 

T5 Brightness temperature of VLBI core component 

R pc Ratio of parsec-scale core to remaining flux density 

Rkpc Ratio of kiloparsec-scale core to remaining flux density 

(3 a pp Fastest measured apparent speed of jet component 

APA Difference in position angle between pc-and kpc-scale jet 



Polarization Properties 
m op t Optical V band linear polarization 

m 5 

TO15 

|xi5 - PA pc 



Mean linear polarization at 4.8 GHz 
Mean linear polarization at 14.5 GHz 
Offset of integrated 14.5 GHz \ w.r.t. pc-scale jet PA 



K 



Deg. 



70 

% 
% 

Deg. 



Tingay et al. (2001) 

This paper 

Various 

Various 

Various 



Impey et al. 1991 
Aller et al. 1992 
Allcr et al. 1992 
Aller et al. 1992 



a http://www. astro. lsa.umich.edu/obs/radiotel/umrao. html 



Note. — Data on properties marked "Various" were gathered from multiple sources (see §3). 



Table 2 



General Properties of Pearson-Readhead Sources 



T A TT 
IAU 

Name 


Other 
Name 

[2] 


Opt. 
Type 

[3] 


Morph. 
[4] 


IDV 

Class 

[5] 




z 

[6] 


EW 

[7] 


Lc 
r< 

l f 


,pc 

'J 


L 5 
[9] 


log 
Rpc 
[10] 


L c ,kpc 
[Hi 


L ex t 
[12] 


log 

Rkpc 

[13] 


Rcf. 
[14] 


L xra y 

[15] 


Ref. 

[16] 


FS-PR sample 


objects 
















































0016+731 




Q 


2 


I 


1 


.781 


17. 





26 


.99 


27 


93 


-0 


.89 


27. 


90 


<25 


.57 


>2 


.33 


1 


20 


.97 


8 


0133+476 


OC 457 


HPQ 


C 


II 





.859 


10 


4 


27 


46 


27. 


53 





.75 


27. 


34 


<24 


.26 


>3 


.08 


1 


20 


.96 


8 


0212+735 




HPQ 


C 





2 


.367 


15. 


.0 


28 


.25 


28. 


78 


-0. 


.37 


28. 


56 


<25 


.59 


>2 


.97 


1 


22 


01 


8 


0316+413 


3C 84 


RG 


2 








.017 


417. 


5 


24 


.07 


25. 


22 


-1. 


.12 


25. 


05 


24. 


58 


0. 


47 


2 


18 


.00 


9 


0454+844 




BL 


C 







.112 


<1. 









25 


.30 






25 


09 


<22 


.71 


>2 


.38 


3 


18 


.01 


10 


0723+679 


3C 179 


LPRQ 


2 







.844 


39. 


.2 














26 


81 


27 


.76 


-0. 


95 


4 








0804+499 


OJ 508 


HPQ 


1 


II 


1 


.432 


13. 


.0 






27 


90 






27 


53 


25 


.75 


1. 


.78 


3 


21 


.23 


8 


0814+425 


OJ 425 


BL 


1 


II 





.245 


1. 


.0 


25 


.83 


26. 


19 


-0 


11 


26. 


38 


25 


.07 


1. 


.30 


3 


18 


.95 


10 


0836+710° 


4C 71.07 


Q 


2 





2 


.180 


17. 


6 


27 


.88 


28. 


51 


-0. 


.52 


28. 


67 


27 


.51 


1. 


16 


3 


22 


.83 


8 


0850+581 


4C 58.17 


LPRQ 


2 




1 


.322 


25. 


9 






27. 


69 






27 


48 


26 


.73 





75 


4 


20 


.79 


11 


0859+470 


4C 47.29 


Q 


2 





1 


.462 


30. 





27 


39 


27. 


95 


-0. 


.42 


27. 


96 


27. 


55 


0. 


.40 


3 


20 


.70 


11 


0906+430 


3C 216 


HPQ 


2 


I 





.670 


3. 


.1 


26 


.74 


27. 


32 


-0. 


.45 


26. 


92 


27. 


22 


-0. 


.30 


5 


20 


.23 


8 


0923+392 


4C 39.25 


LPRQ 


2 








.699 


58. 


3 


26 


.52 


28 


11 


-1. 


.58 


27. 


52 


27 


.09 





.43 


3 


20 


.99 


8 


0945+408 


4C 40.24 


LPRQ 


1 


I 


1 


.252 


15. 


9 






27 


.78 






27 


69 


26 


.86 
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.83 
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20 


.93 


8 


0954+556° 
0954+658° 


4C 55.17 


HPQ 
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.900 


7. 


.1 
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.61 
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BL 
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.23 
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10 
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4C 41.32 
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.550 
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6 


27. 
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4C 38.41 


HPQ 
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.807 
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.58 


28. 


36 
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.70 
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27 


26. 
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1. 


47 


3 
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.53 
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BL 
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.00 
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10 


1928+738 


4C 73.18 


LPRQ 
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.302 


398. 
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.98 
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25 
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OV 591 
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OW 637 


RG 
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BL Lac 
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VSOP-PR sample objects 
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09 


27 


73 


28 
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42 
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08 


9 



a Mcmbcr of third EGRET catalog (Hartman ct al. 1999) 



Note. — All quantities are calculated assuming h — 0.65, q — 0.1, and A — 0. Luminosities arc given in W Hz" 1 . Columns are as follows: (1) IAU 
source name; (2) Alternate name; (3) Optical classification (see §3.1); (4) Kiloparsec-scale radio morphology, where C — compact, 1 — one-sided, 2 — 
two-sided, and H — halo; (5) IDV classification (see §3.2); (6) Rcdshift; (7) Source frame equivalent width of widest permitted line [A]; (8) Luminosity 
of parsec-scale core component at 5 GHz; (9) Total (single-dish) luminosity of source at 5 GHz; (10) Ratio of parsec-scale core flux density to remaining 
(non-core) flux density at 5 GHz in source rest frame; (11) Luminosity of kiloparsec-scale core component at 1.4 GHz; (12) Luminosity of kiloparsec-scale 
extended structure at 1.4 GHz; (13) Ratio of kiloparsec-scale core to extended flux density at 1.4 GHz, in source rest frame; (14) Reference for columns 
4, 11 and 12; (15) X-ray luminosity at 0.1 keV; (16) Reference for X-ray data. 



References. — (1) Xu et al. 1995; (2) Pedlar ct al. 1990; (3) Murphy, Browne, & Perlcy 1993; (4) Reid et al. 1995; (5) Antonucci & Ulvestad 1985; (6) 
Cassaro et al. 1999; (7) Browne &c Pcrley 1986; (8) Sambruna 1997; (9) Hardcastle & Worrall 1999; (10) Urry ct al. 1996; (11) Perlman et al. 1998; (12) 
Comastri et al. 1997; (13) Neff, Roberts, & Hutchings 1995; (14) O'Dea, Barvainis, & Challis 1988. 
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Table 3 



Jet Properties of Pearson-Readhead Sources 



Source 


Pa 




PA pc 


PA kpa 


APA 






Name 


(vfc) 


Ref. 


(deg) 


(degT 


(deg) 


Ref. 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


r o-r^rt sample 


objects 














0016+731 


17.3 


2 


132 


169 


37 


22. 


25 


0133+476 






330 






22. 




0212+735 


8.6 


4 


121 






22, 




0316+413 


0.08 


5 


176 


162 


14 


22, 


26 


0454+844 


1.1 


6 


177 






22, 




0723+679 


8.6 


8 


256 


266 


10 


22. 


38 


0804+499 






127 


201 


74 


22, 


27 


0814+425 






103 


50 


53 


22 


27 


0836+710 


21.7 


9 


201 


198 


3 


22' 


29 


0850+581 


7.6 


10 


227 


151 


76 


22, 


30 


0859+470 






357 


335 


22 


22, 


28 


0906+430 


5.3 


H 


151 


244 


93 


22, 


31 


0923+392 


2.3 


12 


93 


76 


17 


21. 


32 


0945+408 






137 


32 


105 


22. 


28 


0954+556 






191 


301 


110 


23, 


28 


0954+658 


12.2 


13 


326 


224 


102 


22. 


33 


1 694+41 6 






261 


351 


90 


1, 


34 


1633+382 


15.9 


14 


279 


176 


103 


21, 


27 


1637+574 






200 






22, 




1641+399 


29.3 


15 


284 


328 


44 


21. 


35 


1642+690 


14.0 


16 


158 


186 


28 


22, 


34 


1652+398 


5.6 


17 


160 


45 


115 


22, 


36 


1739+522 






204 


260 


56 


22, 


28 


1749+701 


10.9 


2 


275 


21 


106 


22. 


25 


1803+784 


3.2 


13 


291 


193 


98 


22. 


27 


1807+698 






252 


261 


9 


22, 


36 


1823+568 


21.2 


14 


196 


182 


14 


22. 


34 


1928+738 


11.5 


4 


166 


190 


24 


21, 


27 


1954+513 






307 


338 


31 


22, 


24 


2021+614 


0.19 


19 


205 






19, 




2200+420 


13.0 


20 


209 






22, 




2351+456 






321 


348 


27 


22, 


39 


VSOP-PR sam 


pie objects 












0153+744 


3.8 


3 


68 






22. 




0711+356 


< 0.81 


7 


329 


294 


35 


22, 


27 


1828+487 


10.6 


18 


311 


311 





22. 


37 



Note. — Columns arc as follows: (1) IAU source name; (2) Maxi- 
mum component speed in units of c, assuming h — 0.65 and q a — 0.1; 
(3) Reference for apparent speed; (4) Parscc-scalc jet position angle; 
(5) Kiloparscc-scalc jet position angle. (6) Difference in parsec- and 
kiloparscc-scalc jet position angles; (7) References for images used in 
columns 4 and 5, respectively. 



References. — (1) Lister et al. 2001; (2) Schalinski et al. 1992; (3) 
Hummel et al. 1997; (4) Witzel et al. 1988; (5) Dhawan, Kcllcrman, 
& Romney 1998; (6) Gabuzda et al. 1994; (7) Conway et al. 1994; 
(8) Porcas 1987; (9) Marchenko-Jorstad et al. 2000; (10) Barthel et al. 
1986; (11) Paragi et al. 2000; (12) Fey, Eubanks, & Kingham 1997; (13) 
Gabuzda & Cawthornc 1996; (14) S. G. Jorstad, private communica- 
tion; (15) Ros, Zensus, & Lobanov 2000; (16) Pearson et al. 1986; (17) 
Giovannini et al. 2000; (18) Polatidis & Wilkinson 1998; (19) Tschager 
et al. 1999; (20) Denn, Mutel, & Marscher 2000; (21) Lister & Smith 
2000; (22) Lister et al., in preparation; (23) A. P. Marscher, private 
communication; (24) Kollgaard, Wardle, & Roberts 1990; (25) Xu et 
al. 1995; (26) Pedlar et al. 1990; (27) Murphy, Browne, & Perley 1993; 
(28) Reid et al. 1995; (29) Hummel et al. 1992; (30) Shone. Porcas, fc 
Zensus 1985; (31) Taylor, Go, & O'Dea 1995; (32) Marscher et al. 1991; 
(33) Stanghcllini ct al. 1997; (34) O'Dea, Barvainis, & Challis 1988; 
(35) Kollgaard, Wardle, & Roberts 1989; (36) Cassaro ct al. 1999; (37) 
Liidke et al. 1998; (38) Akujor 1992; (39) Neff, Roberts, & Hutchings 
1995. 
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Table 4 

Correlations Among Observed Properties of the FS-PR Sample 



Property 


Property 


N 


r 




P 




Tz 


Pz 




Notes 


(!) 


(2) 


(3) 


(4) 




(5) 




(6) 


(7) 




(8) 


r 

^ext 


T> 

*^kpc 


30 


-0.56 


9.7 


X 


1 fi- 
ll) 


-4 


-0.60 


4.9 x 10" 


-12 


Expected 


L^pc 


EW 


24 


-0.49 


8.1 


X 


1 fi- 
ll) 


-2 


-0.50 


1.8 x 10" 


-7 


Impey et al. (1991) 


L^ 


■Lc,kpc 


28 


0.91 


1.1 


X 


1 fl- 

1U 


-9 


0.77 


7.6 x 10" 


-7 


Akritas & Siebert (1996) 


v lo 


Vk 

Y 


31 


0.62 


9.6 


X 


10" 


-5 


0.63 


7.2 x 10" 


-5 


Aller et al. (1992) 


■Lext 


^320 


28 


0.48 


2.4 


X 


10- 


-2 


0.49 


9.2 x 10" 


-4 


Expected 


Wlopt 


EW 


30 


-0.42 


1.2 


X 


10- 


-1 


-0.40 


2.1 x 10" 


-2 


Impey et al. (1991) 


L 5 


L X ray 


25 


0.79 


3.6 


X 


10- 


-6 


0.52 


2.6 x 10" 


-2 


Browne & Murphy (1987) 


L pt 


Lxray 


26 


0.67 


1.4 


X 


10" 


-4 


0.52 


2.6 x 10" 


-2 


Fabbiano et al. (1984) 


Rpc 


m op t 


23 


0.50 


8.2 


X 


10" 


-2 


0.52 


4.5 x 10" 


-2 


Impey et al. (1991) 


APA 


EW 


26 


-0.39 


5.5 


X 


10- 


-1 


-0.40 


2.0 x 10" 


-1 


Appl et al. (1996) 


L 5 


L320 


27 


0.69 


5.1 


X 


10- 


_5 


0.42 


2.2 x 10" 


_ y 


Zirbel & Baum (1995) 


L 5 


Lc,pc 


24 


0.75 


2.5 


X 


10- 


_5 


0.41 


4.8 x 10" 






Lc,pc 


Lc.kpc 


23 


0.76 


4.6 


X 


10- 


_5 


0.42 


4.9 x 10" 


_ y 


Expected 


mi 5 


m 5 


01 




7 9 


X 


10- 


-1 




71 v in- 
/ . 1 x iu 


-1 




L e xt 


L 5 


28 


0.33 


1.1 








0.25 


9.3 x 10- 


-1 


Expected 


Lc,kpc 


-Lxray 


25 


0.71 


2.5 


X 


10- 


-5 


0.41 


1.2 




Browne & Murphy (1987) 


Lc,pc 


-Lxray 


20 


0.72 


1.0 


X 


10- 


-3 


0.39 


1.6 




Bloom et al. (1999) 


L op t 


Lc.kpc 


30 


0.60 


3.2 


X 


10- 


-4 


0.31 


1.7 




Feretti et al. (1984) 


v 5 


«5/15 


31 


0.30 


1.6 








0.30 


1.7 




All j_ 1 /1 AArt\ 

Aller et al. (1992) 


L c ,kpc 


£320 


28 


0.63 


2.3 


X 


10- 


-4 


0.31 


1.9 




Fabbiano et al. (1984) 


m pt 




29 


0.33 


1.3 








0.31 


2.0 




Valtaoja et al. (1992) 


Correlations likely due to 


rcdshift 


bias 


















L pt 


ftapp 


20 


0.38 


2.0 








0.27 


4.4 






L op t 


■Lc,pc 


24 


0.58 


8.4 


X 


10- 


-3 


0.28 


6.0 






L pt 


L 5 


30 


0.54 


2.4 


X 


10- 


-3 


0.20 


13 






Lc,pc 


L320 


22 


0.56 


2.8 


X 


10- 


-2 


0.19 


22 






L320 


Lxray 


23 


0.48 


1.4 


X 


10" 


-1 


0.07 


65 






L pt 


L32O 


29 


0.36 


5.4 


X 


10" 


-1 


0.02 


86 






Rejected correlations 






















Ixi5 -PA pc \ 


EW 


20 


0.40 


1.5 








0.42 


3.8 x 10" 


-1 


Due to single outlier 


m op t 


"5/15 


29 


0.32 


1.6 








0.32 


1.6 




Due to single outlier 


m c 


Papp 


19 


0.39 


1.9 








0.38 


1.8 




Due to single outlier 



Note. — Columns are as follows: (1) and (2) Source property (see Table 1); (3) Number of data pairs; 
(4) Kendall's tau coefficient; (5) Probability of correlation arising by chance in per cent; (6) Kendall's 
tau coefficient, with redshift partialed out; (7) Probability of correlation arising by chance, with redshift 
partiallcd out; (8) Comments and/or reference for previous detection of correlation, if applicable. 
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Table 5 

Two-Sample Tests for Various Subclasses of the FS-PR Sample 





QSO 


HPQ 


IDV II 


IDV I & II 


IDV II 


2-sided 




vs. 


vs. 


vs. 


vs. 


vs. 


vs. 


Property 


BL 


BL 


IDV & I 


IDV 


IDV 


Others 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


z 


0.070 


1.9 










EW 


0.0053 


0.017 




0.13 


0.16 


0.021 


«5/15 








0.46 


0.40 




v 5 






0.17 




0.077 




v 15 














L320 


0.16 


0.50 




0.90 


1.2 




L 5 


0.035 


0.48 




0.25 


1.7 




L pt 


1.4 












Lxray 


0.31 












Lc,pc 














Lc,kpc 


0.078 


0.098 










Lext 


1.72 










0.011 


Rpc 






0.19 




0.082 


0.46 


Pkpc 












0.098 


m op t 


0.78 




0.35 


0.27 


0.22 




APA 












0.54 


T b 








1.98 







ftapp 

m 5 
mi 5 

|xi5 - PA pc \ 



Note. - All values represent the probability (in per cent) that the sub- 
samples were drawn from the same population, based on the source property 
in question. Those values not listed exceed 2%. Columns are as follows: (1) 
Source property (see Table 1); (2) Quasars vs. BL Lacs; (3) High-optical po- 
larization quasars vs. BL Lacs; (4) IDV class II AGNs vs. the combined IDV 
classes & I; (5) Combined IDV class I & II vs. class 0; (6) IDV class II vs. 
class 0; (7) Sources with two-sided kiloparsec-scale morphology vs. all others. 



